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1. INTRODUCTION 

The Calogero-Sutherland-Moser model (CSM) [1-3] has attracted much 
attention in recent years due to its relation to a wide variety of interest- 
ing problems. Some examples are random matrix theory 0], quantum spin 
chains with long-range interactions [^, generalized exclusion statistics [6-12], 



Gaussian conformal field theories [13|, edge states in a quantum Hall system 



IJ], and nonlinear internal waves in a stratified fiuid |]T5 . 

The CSM has been known to be exactly solvable and integrable, both 
classically and quantum mechanically, for quite some time [1-3, 16]. However 
detailed investigations into its collective properties have begun only recently 
[17-19]. A collective field theory to study the excitations of a superfiuid, 
as well as the ground state of a condensed Bose-Einstein gas was developed 
long back [^. In such a theory, the fundamental coordinate is the density 
field [21-24]. For the CSM, the results obtained so far include the ground 
state energy of the model placed in a harmonic oscillator potential, waves of 
arbitrary amplitude for strong coupling, and isolated solitons on an uniform 
background density. 

In our paper, we will study essentially the same aspects but in more detail 
and for arbitrary coupling, thereby generalizing the earlier results in several 
ways. Wherever appropriate, we will compare our results with those obtained 



earlier by other methods |]TI|, ^. This will illustrate that certain properties 
of the model can be derived more easily and generally from collective field 
theory. These properties include the dispersion relation for small amplitude 
and long wavelength sound modes, the low temperature specific heat, and the 



two-point correlation function. The collective field theory yields a dispersion 
relation for the sound mode that terminates exactly in the second order of 
the wave number, which is adequate only for small wave numbers. This gives 
the correct nonos dilatory behavior of the two-point correlation function for 
asymptotically large distances, and it fails at shorter distances. The collective 
field theory formulation is thus seen to be a useful description for the low- 
energy, or small wave number excitations of the CSM system. 

2. CALOGERO-SUTHERLAND MODEL AND COLLECTIVE 

FIELD THEORY 

The simplest form of the CSM consists of particles on a line which interact 
pairwise through an inverse- square potential. The model can also be defined 
on a circle with periodic boundary condition |0]; the two versions of the 
model have identical physical properties in the thermodynamic limit in which 
the the number of particles N and the length L of the line (or circle) are 
simultaneously taken to infinity keeping the particle density po = N/ L fixed. 
The Hamiltonian for particles on a line is given by 

^ 2m m ^- ixi- x,)2 ' ^ ' 

where the dimensionless coupling A > 0. To make the problem well-defined 
quantum mechanically, we have to add the condition that the wave functions 
\E' goes to zero as |xj — Xj\^ whenever two particles i and j approach each 
other. For A = and 1, the model describes free bosons and free fermions 
respectively. Since the two-body potential is singular enough to prevent par- 
ticles from crossing each other, we can choose the wave functions to be either 



symmetric (bosonic) or antisymmetric (fermionic). The energy spectrum is 
the same in the two descriptions. 

Let us briefly summarize some of the exactly known results for this model 
[0, |T^. If Eq denotes the ground state energy, then the chemical potential 
at zero temperature is given by /x = OEq/ON in the thermodynamic limit. 
This takes the form 

In a fermionic description of the model, it is natural to define a Fermi mo- 
mentum 

Po = irhpo . (3) 

(We should point out that some papers in this field find it more convenient 
to define the Fermi momentum to be vr/iApo)- The low-energy excitations of 
(IID are known in detail. They can be thought of as being made up of particle 
and hole excitations. Let us first define the sound velocity by the relation 



Then 

nhXpQ 



(5) 



m 
It is known that the particle excitations necessarily have \p\ > Pq, with the 

dispersion 

ep(p) = ^ (bl - Po) { \P\ + (2A - l)po ] + /i • (6) 

The hole excitations have |p| < po, with the dispersion 

e^(p) = ^(^0-/) - /^- (7) 



If we define group velocities v = hdu/dp, we find that particles have \v\ > Vs, 
while holes have \v\ < Vg- The sound velocity in (j^) can be obtained by 
considering a sound mode to be made up of a particle with energy-momentum 
(ei,pi) and a hole with energy-momentum (€2,^2)- Then Vg = (ei + e2)/(pi — 
P2) in the limit pi,P2 — > Po- 

The collective field theory for the CSM is obtained by changing variables 
from the particle coordinates Xi to the density field p{x) defined as 



As emphasized in reference [^ , such a change of variables is meaningful only 
if the particle number N —>■ 00. We therefore have to check at various stages 
whether the results obtained from collective field theory for finite values 
of A^ are indeed correct for the model defined in (|l|). For this reason, we 
will compare the collective field theory results with those obtained by other 
methods whenever possible. 



After changing variables, the quantum Hamiltonian takes the form ||22|, ^ 



H = — / dx 

2m 



2 , _2\2 



4 p J 

where d = d/dx, and pn is proportional to the Hilbert transform of p defined 
as the principal part integral p7| . 



Ph{x) = dy p{y) , 

J X — y 

P 1,1 1 , 

where = - lim ( + ) . (10) 

X — y 2 e^o X — y + ie x — y — it 



The field W is canonically conjugate to p, and they satisfy the equal-time 
commutation relation 

[ p{x) , h9{y) ] = in 5{x - y) . (11) 

We may therefore set 9{x) = —i6/6p{x), and try to find eigenstates of the 
Hamiltonian (|^) as functionals of p(x). Although this is sometimes possible 



2^, it is generally very difficult to find exact eigenstates. 

We will therefore take the simpler route of studying (^ classically. For 
this purpose, it is useful to rewrite the collective theory in terms of a fluid. 
Let us introduce a complex field 

tlj{x,t) = ^p{x, t) e*^(^'*) . (12) 

We define a Lagrangian L = J dxC, where the Lagrangian density 

£ = y (^> - ip*ip) - — dip*dip - U[p{x)] 

where UW = i^l f i!^ / + A(A - 1) ,„gp + ^^^ '^"'^ 



P 



(13) 



and a dot denotes d/dt. From (|l^) we see that p and h6 are canonically 
conjugate to each other, and we can recover the Hamiltonian (^ by the 
usual methods. 

It is interesting to note that the Lagrangian (p!3|) is quadratic in ip and 
ip*, and is therefore noninteracting for A = (free bosons). This is under- 
standable because the collective field theory is a bosonic theory, as is clear 

6 



from (|TT]). On the other hand, the collective theory is interacting for A = 1 
(which is free in terms of a fermionic theory). 

We now proceed to study the theory ([T3|) classically. At a classical level, 
Eq. d^) (for the static case) may be regarded as the energy density functional 
of the single-particle density p{x). This is analogous to the density functional 
theory of a correlated many-particle system, a highly sucessful formalism in 
many branches of physics [^. Since we wish to study the system in an 
external potential V{x) in Section 3, let us add 

f dx[fx- V{x) ] p{x) - pN (14) 

to the Lagrangian, where fi is the chemical potential. The Euler-Lagrange 
equations of motion obtained by extremizing the action S = J dtL are given 

by 

Ti^n^x^ 2 _ ^'A(A - 1) !^!(Az^ r r^f _ 2^ 1 

+ he + — [oef - fi + v{x) = , (15) 

^ lit 



and 



p + — dipde) = . (16) 

m 



In addition, the density must satisfy the constraint 

fdxp = N . (17) 

Eq. (|l^) will be recognized as the equation of continuity since {Ti/m)d6 is 
the velocity field; this can be seen from the expression for momentum given 
below. 

7 



Our system has three conserved quantities, namely, the particle number 
A^, the momentum (if there is no external potential) 

P = - y f dx iij^dip - dip^ip) = h fdxpde, (18) 

and the energy 

^' r , r ,.„.. vr^A^ , , ,,, ^,^ ^^ , (A-1)^ (dpf 



E= ^ Jdx[ p{def + ^p3 + A(A - i)pHdp 



4 p 

dx V{x) p . (19) 

There are probably an infinite number of conserved quantities in addition to 
the three above since our original system is integrable; however explicit field 
theoretic expressions for these other quantities are not known. 

In the absence of an external potential V{x), Eqs. (|I5|) and (|1^) are 



invariant under scaling and Galilean transformations. Under scaling by a 
factor a, we have 





p{x 


t) ^ a p{ax, a'^t) 




e{x 


t) -^ e{ax,aH) , 
p —>■ a^p . 


Under a Galilean transformation by velocity w , 


p(x,t) 


-^ 


p{x-vt,t) , 


e{x,t) 


- 


9{x-vt,t) + — ( 


/^ 


^ 


p . 


Thus 






F 


) _ 


-^ P + mNv , 



(20) 



1 
2 



(21) 



E -^ E + Pv + -mNv^ . (22) 

It is quite remarkable that if ^ = 0, all values of A > 1 are equivalent to 
each other according to Eqs. (|T3[) and (0). Namely, if we redefine 

X = X , i = t{\ — 1) , 
Ap ~ 9 

P = -^ 7 ' ^ 



A-1 ' A-1 ' 

^ - (X^. (23) 

then A — 1 can be completely scaled out of (|15D and ([T6|) . Similarly, all values 
of A < 1 (but not equal to 0) are equivalent to each other; we can carry 
out the same redefinitions as in (pSj), followed hj p ^ —p in order to keep 
p positive. If we redefine the energy functional as ii^ = XE/\X — Ip, we see 
from (H) that 

^2 . . „2 ,^.^2 



E 

2m 



jdx [ p{d~ef + ^p3 ± Pndp + ^ ] , (24) 



where the ± signs are for A > 1 and A < 1 respectively. Thus it is sufficient 
to study the collective field theory for just two values of A, one less than 
1 and the other greater than 1. This property of the collective field theory 
clearly shows that it is a rather coarse description of the CSM. This is to 
be contrasted with the exact solution of the model (0) some of whose fea- 
tures (for instance, the dynamical correlation functions [^) are sensitively 
dependent on number theoretic properties of A. 

In passing, it may be noted that formally the scaled energy density given 
by Eq. ( p4D is of the same form as the so-called Madelung fiuid [^ , which 



is a hydrodynamical description of the one-particle Schrodinger equation. In 

9 



this picture, the first term is the classical kinetic energy of the fluid, the 
next two represent the potential energy, and the last term arises from the 
quantum kinetic energy. The latter gives rise to the Bohm potential [^ 



m 



the equations of motion. This interpretation also holds if there are A^ particles 



in the same quantum state. The normalization of p in Eq. (]2^) , however, is 
not N . To pursue this line of thought more carefully, it is necessary to modify 
transformations (p3D, and scale the x-coordinate to demand J dxp = N. The 
Bohm term in ( p4|) remains unaffected, but the interaction terms become 
A-dependent. We will not elaborate further along these lines. 

We will now study various solutions of the equations of motion (|T^ and 
(|16D, with an external potential in Section 3 and without an external poten- 
tial in Sections 4-6. We are interested in two kinds of solutions, (a) static 
solutions in which p depends only on x and 6 = {in particular, the ground 
state is always of this form), and (b) time-dependent solutions in which p 
and 6 depend on both x and t. 

3. GROUND STATE IN AN EXTERNAL POTENTIAL 

For any external potential, Eq. (|r^) gives the exact quantum ground 
state energy and density if A = 0. In that case, let "^olx) and cq denote 
the exact one-particle ground state wave function (normalized to unity) and 
energy obtained by solving the Schrodinger equation with a potential V{x). 



Then the solution of ([T5|) and (|T^) is given by 

p{x) = N\^o{x)\\ 
Eo = Neo . (25) 

10 



The question therefore is how well collective field theory does for nonzero 
values of A. 

To begin with, let us consider the case of a simple harmonic potential, 
with 

V{x) = -muj^x'^ . (26) 

This turns out to be a rather special case because the ground state energy of 
the collective field theory can be found exactly. Since this is a static solution 
with 6' = 0, we can use the principal part identity 

P P P P P P 

— ^— + — ^— + — ^— = TiH{x-y)5{x-z) , (27) 

X — y X — z y — z y — X z — x z — y 

to write (|TU|) as a perfect square 



E = — J dxp{Xpn + — --—^) +^[XN' + {1-X)N]. (28) 



h f j^ I' ^ \ — 1 dp muj ^^2 ^ huo 

P 
Thus if p satisfies 



X — 1 dp mu ,„„. 

Ap„ + ^^ = ^.. (29) 

then it minimizes (^) and is therefore a solution of the equation of motion 
(|T5|). Further, the ground state energy follows from (pHD, 

Eo = ^ [\N' + (1 - \)N] . (30) 

This is in fact the exact answer for the Hamiltonian (|1|). 

Eq. (^Of ) for the density can be solved analytically only if A = or 1. We 

get 

p = N ( ) exp ( — mujx'^ Ih ) if A = 

mw / 2Nh , \i/2 ^ ^ ^ ^ 

= -— ( xM if A = 1 . (31) 

Tin ^ mu ' 

11 



We can show analytically that the collective field theory density has a Gaus- 
sian tail of the form 

if A < 1, and has a sharp cutoff x = ±Xq beyond which p vanishes if A > 1. 
For large values of A^, we can also show that the second term on the left 
hand side of (^) is generally much smaller than the first term; the scaling 
argument is indicated below. If we ignore the second term altogether, we get 
the leading behavior of p to be a semicircle for all nonzero values of A, 

mw / 9 9 \i/2 



P 



Here Xq is defined by 



Ixq — X j for \x\ < Xo, 



TihX 

for |x| > Xq . (33) 



Xo = ( ) • 34 

The relations (|33|) and (p4D are identical to the Thomas-Fermi result obtained 
in [|T^, and Xq is just the classical turning point. Note that the form of 
p in (^) is essentially a statement of exlusion statistics for the GSM; the 
occupation number in each state in phase space dxdj) = 2Tth is given by 1/A. 
One should, however, be wary of using the expression (^) for the density 
p{x). For example, if we indiscriminantly substitute this p{x) in the static 
energy density functional 



dxV{x) p + — dx[^—p^+ \{\-l)pHdp + 



P 
(35) 
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the integrals with V{x) and p^ on the right-hand side together yield the 
correct A^^-dependent term in Eq, and the third integral gives the right A^- 
dependent term (see Eq. (pOD, but the last integral involving {dpY / p diverges. 
It may be easily checked, however, that this divergent term goes like A^°, i.e., 
of order 1 in the large-N expansion. Such terms in the expansion will be 
dropped. 

We should point out that the 1/A^ expansion of p within the collective 
field theory cannot be taken too seriously; we recall the cautionary remarks 
following Eq. @. For instance, the absence of a Gaussian tail if A > 1 is an 
artifact of collective field theory. If A^ is finite, the ground state of (|I]) has a 
Gaussian tail for all values of A; this can be seen from the exact expression 

^0 [x^] ~ n l^i -^il^ exp [ - — - ^ xH • (36) 

i<j ^^ i 

On the other hand, we can generally trust the next to leading term in the 

energy given by collective field theory. We have already seen this for the 

harmonic oscillator potential, where the N"^ and A^-dependent terms both 

came out correctly. We will now show this for a somewhat more general 

class of potentials. 

We formally define the first two terms in a 1/A^ expansion as follows. We 

assume that the ground state energy and chemical potential have expansions 

of the form 

p = /i(°) + /i(^) , 
Eo = 4°) + 4') . (37) 

where Eq /Eq and p^-^^ / p^^^ are of order 1/A^. For the density, we have an 

13 



expansion of the form 

p{x) = p^''\x) + p^'-\x) . (38) 

We will state the procedure for obtaining the expansion (|38| ) shortly. But we 
can point out immediately that it is not a 1/A^ expansion for all positions x; 
although p^^\x) / p^'^'> {x) will generally he order 1/A^, that will not be true near 
the turning point xq. The first term p*^"^ in ( p8| ) is defined by considering only 
the p^ and Vp terms in the energy functional (^) and hence in Eq. ([T5|) . 
The leading term in the chemical potential /i^^-* and the turning point xq 
are then fixed by the particle number constraint (|T^. Next, we define the 
term p^^-* in (|38|) by considering the p^, Vp and pndp terms in (|35| ) and the 
corresponding terms in ([T5|). Once again, p^^^ is fixed by the constraint ([T7|). 
We will not go beyond the two leading terms in Eqs. ( pTf ) and (pS]), and will 
therefore not need to consider the {dpY / p term in (p5|); this term actually 
diverges even more severely for p^^^ than for p^^\ 

As a specific example of the 1/A^ expansion, let us now consider a con- 
fining potential of the power-law form 

^'W - i^ ( T )' ' (-' 

where p > 0, and a is a measure of the width of the potential. (A harmonic 
potential corresponds to the case p = 2). For such a potential, we can 
prove that the Vp + p^, pndp, and {dpY / p terms in the energy functional 
(|35D are successively of higher order in 1/N . To show this, let us define the 
dimensionless variables 

1 X 

X 



Ar2/(P+2) a 
14 



p{x) 



1 



a p{x) 



so that J dxp = 1. Then Eq. (|35|) takes the form 



(40) 



E 



2ma'^ 



^(3p+2)/(p+2) /^~ 



L^IP ;^ 



xl'^ p + 



+ 



TT^A^ 



. A(A-l). .^ 
P + .r Pndp 



3 ' ' N 

(A - 1)2 {dpr 



(41) 



4A^2 p 

This justifies the form of the 1/A^ expansion given in the previous paragraph. 
Following that procedure, the leading order terms in p and p are found to be 



p 



(0) 



Xo\P 



p 



(0) 



2ma'^ ^ a 

-^ [ 2m(/i(o) - V{x) ) ]^/2 if |x| < xo 



(42) 



From the constraint (|T7|) we find 

ttAA^ 



'Xo\l+p/2 

^ a ' 

h 



2/i 
1 







dy Jl - yP 



r(i + ^)r(|: 

ni + i) 



(43) 



The leading terms p and Vp in the energy ( pSf ) then give 



E, 



(0) 



xo 



dx [ 2m(/i(°) - V{x) f^ (^(°) + 2V{x)) 

jY{3p+2)/{p+2) _ /^^\ 



Svr/iA Jo 

h'^ p + 2 /7rA\2p/(p+2) 
2ma2 3p + 2 ^2/^/ 

We now go to next order in 1/A^ by including the terms in (|15|) and (p5|) 
which contain the Hilbert transform pn] we get 



n^h^X^ 



m 



p(o)^(i) 



^2A(A - 1) 



(0) 



m 



dp'n 



p 



(1) 



(45) 
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The structure of ( ^Sj) shows that p^^^ must be taken to be zero outside the 
turning point Xq, just hke p^^\ Since p*-*'-' is normahzed to A^, we fix p^^^ using 
the constraint 

/ dx p(^) = , (46) 

J—XO 

and we then determine Eq by including the pndp term in (|35|). Interestingly, 
we find that p*^^) = for all A for the simple harmonic case (p = 2); that is 
why we get the energy correct to order N by just substituting p^^' in (|35|) . 
However p*^^-* is not zero for a general value of p. We will omit the final 
expression for Eq for general A, and will now specialize to A = 1 where we 
can compare with the results of a WKB approximation. (Although this is a 
free fermion theory, it is an interacting bosonic theory. Hence agreement at 
A = 1 is a nontrivial check of collective field theory). We find that both p^^^ 
and p^^-'(x) are zero for A = 1. Hence there is no correction to Eq at the next 
order after (|4^), i.e., at order N'^p/(p+'^\ We will now show that this result 
agrees with WKB. 

If Cn denote the single-particle energy levels obtained by solving the 
Schrodinger equation in the potential V{x), then the exact ground state 
energy at A = 1 is given by 



N-l 



Eo = E en . (47) 

n=0 

For large n, we can obtain the two leading order terms in e„ using the WKB 

formula 

dx[2m{en- V{x) ) f^^ = ( n + - ) nh , (48) 

-x„ ^ 2 ^ 

where Xn denotes the classical turning point for energy e„. We thus obtain 
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the expansion 

h'^ / 7m \ 2p/{p+2) p 



2ma2 V2Jiy ^ ^ (p + 2)n ^ ' ' ' ^ ' ^^^^ 

On substituting this in (0), we find that -Eq is indeed given by (|^) and that 
there is no correction to the next order in 1/A^. 

4. SMALL AMPLITUDE WAVES, CORRELATION FUNCTIONS, AND 

SPECIFIC HEAT 

In this Section, we will study the small amplitude density fiuctuations 
about an uniform background density po- We will show that these exhaust 
the low-energy excitations upto some order, both in a sum rule and in the 
low-temperature specific heat. 

For an uniform density po, the chemical potential is given by (0) or (p!5|) 
to be 

„2t2x2„2 

^--^- (50) 

Let us now study ([T5|) and (|16]) to linear order in an amplitude a << 1. We 
assume 

P = Po + ctPo cos(A;x — ut) , 

muj 
9 = a — — sin(A;x — cot) , (51) 

ilrh 

where k denotes the wave number; the second equation in (|5l|) follows from 



the first due to the equation of continuity ([iq ). Eq. (^) then yields the 
dispersion relation 

Tr^Apol^l (A - l)hk'^ 

17 



The sound velocity is given by the group velocity duj/dk at A; = 0; the result 
agrees with the exact value given in (j^). Note that (^) gives the correct 
single-particle dispersion for A = 0, as expected for a free boson theory. 
We see that the dispersion (|52D vanishes not only at /c = 0, but also at 



1*1 = f^ (53) 

if A > 1. However the latter point where uj vanishes seems to be an artifact 
of collective field theory; it does not agree with known results. For instance, 



reference [25| defines a dispersion relation called the "Feynman spectrum" as 
follows. Consider the dynamical correlation function and its Fourier trans- 
form 

G{x,t) = (p(x,t)p(0,0)) - pI, 

S{k,ij) = -^ [ dx [ dt G{x,t) e-'^^'-'-'^'K (54) 

zvrpo J J 

S{k,u!) can be represented in terms of all the states of the system \n) with 

energies En as 

-' ' n 
N 

Pk = Y: e~''^- . (55) 

n=l 

One can then define various moments of S{k, u) as 

J„(fc) = /dcjcj" S{k,uj) , (56) 

where Ii{k) = k"^ /2m. The Feynman dispersion is defined as 

.,(*) - II . (57) 



Now this dispersion is known from random matrix theory [^ ^ for three 



special values of A = 1/2, 1 and 2. For these three values, it is found p5 
that ujp{k) agrees with (|52D upto second order for k close to 0. We therefore 
believe that ujp{k) is given by ( |5^ ) to order k'^ for all values of A. However, 
the agreement between ujp{k) and (^) does not persist to higher orders in 
k even for the three special values of A; in particular, upik) does not vanish 
at any nonzero values of k, although it does have a roton-like minimum at 
\k\ = 271 po for A > 1 1^. This discrepancy between upik) and ( ^2|) seems 



to indicate that collective field theory cannot be trusted for large values of 
the wavenumber; it seems to work only upto order k"^. 

Our statement that the low-energy dispersion is correctly given by (^) 
upto k"^ near k = also agrees with the known low-temperature specific heat 



of the GSM to second order in the temperature T [^ . We can compute the 
free energy per unit length / from (|52|) taking the sound modes to have zero 
chemical potential. Thus 

^ In ( 1 - e-^'^-'' ) , (58) 

-oo ZTT 

where f3 = l/ksT. After evaluating this, we can obtain the specific heat per 

unit length Cy = —Td^f/dT'^ to second order in T. We find 

Cy = !^ + 6C(3)(A-1) klT^ ^ 
3hvs TT mhvl 

oo -| 

C(3) = Y.-,- (59) 

n=l "■ 



This agrees with the result in reference [^ 



We note that the linear terms in (|52D and (|59|) are typical of a system 



whose low-energy and long- wavelength excitations are governed by a confor- 
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mal field theory |T^. The quadratic terms in those two equations indicate 
deviations from conformal field theory which start appearing when the wave- 
length is no longer much longer than the typical particle spacing 1/po- 

Finally, it may be useful to see what we get if we compute the correlation 
function G{x, t) defined in (|5^) by quantizing the collective field theory. Using 
the commutation relation (|l^) and the equations of motion (|15]) and (|l^) 
to linear order, we find that p and 6 have the following second quantized 
expressions, 

'•°° dk 



Po + / — /fc 
i f^ dk 1 



2 J-oo 2n fk 
hpok"^ 
2mujh 



ai. e^e^^-^fc*) -\- a} g-*('=^-'^fe*) 



fljL e^'^'^^~'^kt) _|_ ^t ^-i{kx-uj^t) 



h = C^)"\ (60) 



where 



[ Ofc , 4, ] = 27r 5{k -k') . (61) 

From this we find that 

hpo f^ dk A;2 ^.(^^_^^^j 



G{x,t) = -^ / gHfe^-^fet) . (62) 

2m J-oo 2tt ujk 

If we now use the collective field theory dispersion (^), the integral will 
diverge at the nonzero values of |A;| where u vanishes. We will therefore 
assume, as stated above, that (^) can only be trusted in the region near 
A; = 0. The asymptotic form of G{x,t) at large values of x ± f^t only gets 



a contribution from that region in the integral (]6^) ; further, only the linear 



term in the dispersion (ET) is required to derive the asymptotic expression. 
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We then find that 
G{x,t) 



(63) 



This agrees with the leading nonoscillating term in the exact expression given 
in reference |^8|]. However the exact expression also has oscillating terms; in 



fact, such a term dominates over ([63| ) if A > 1. The fact that collective 
field theory is unable to reproduce these oscillating terms clearly shows its 
limitation. 



5. LARGE AMPLITUDE WAVES 



Following a method given in reference [T^, we can find exact solutions 
which describe waves with arbitrary amplitude. We will consider the cases 
A > 1 and < A < 1 separately. 

For A > 1, the solutions are given by 

(A — 1)A; r sinha 



p{x,t) 



27rA 



mv 
{x,t) = — 



(c+1) 



cosh a — cos k{x — vt) 

cosh a — cos k{y — vt) 



dy 



+ X - - vt 



oo ccosha + sinha — ccos/c(|/ — ft) 

(64) 



We choose the wavenumber k and the parameter a to be positive. The phase 
velocity v in ( |64D satisfies 



where 



2 (A - ifn^^c^ ,2 

^ = \ o/ , 1N2 c' + 2ccotha + l 
4m2(c+l)2 ^ 

1 — cosh a 
c > 



sinha 
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(65) 



(66) 



The average density for this solution is found to be 

We can define the dimensionless amphtude a of the wave to be the fractional 
difference between the maximum and minimum densities. Thus 

Pmax Pmin ^ fnr,\ 

a = = — r-; . (68) 

Pmax + Pmin C Smh a + COSh « 

For A < 1, the solutions are again given by Eqs. ( |64D and (]65|), but 

— 1 — cosh a . . 

c < — . 69 

smh a 



The average density for this solution is the same as in Eq. (|67D , while the 

amplitude a is 

a = — . (70) 

c smh a + cosh a 

The solutions above are characterized by three independent parameters 

which may be considered to be the average density po? the wavenumber k, 

and the amplitude a. If we hold po and k fixed and let a — ^ oo, we recover 

the small amplitude waves discussed in the previous Section. Let us now 

look at the conditions under which the frequency uj = \vk\ can vanish; this 

corresponds to stationary waves. We can see from Eq. (|65|) that u vanishes 

if fc = or if c = 0; the latter is allowed if A > 1 in which case k satisfies 

(|53D. These two conditions for u = are therefore the same as those found 

for the sound modes in Section 4. However we now see that uj also vanishes 

if 

1 — cosh a . . 



sinha 
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for A > 1, or if 

— 1 — cosh a ,„„, 

c = — 72 

smno; 

for A < 1. In these two cases, we get staionary waves with the largest possible 
amplitude a = 1 since pmin = 0. 

The interpretation of all these large amplitude waves, including the new 
kinds of stationary waves in ( fTTj ) and ([72D, in terms of the exact solutions 
of the CSM model remains an open question. It is possible that some of 
the solutions obtained here are peculiar to the collective field theory and do 
not correspond to anything in the CSM. For instance, there are no exact 
solutions of the CSM which have arbitrary nonzero values of k with u = 0. 

In concluding, we would like to mention that the density waves studied 



in this Section and in the previous Section were known earlier [113] for large 



A. In addition, the stationary waves ([7^) for A < 1 were found in reference 



rn]. Our own results describe both stationary and moving waves, and are 



valid for all values of A. 

6. SINGLE SOLITONS 

We will now describe the single soliton solutions of the collective field 
theory |[l^, |19[. Starting from the large amplitude waves in Section 5, we 



can find these solutions for any A different from and 1 as follows. In Eq. 
(|64D , we take the limit k,a —* keeping a/k = b fixed. Simultaneously, we 
let c -^ oo for A > 1 or — cxd for A < 1, keeping po fixed according to (|67|). 
Since the wavelength 2TT/k -^ cxd, we obtain a solution describing an isolated 
lump. We find the following expressions for p and 9 in terms of the width h 
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and velocity v of the soliton. 



p(x,t) 
where v 

T) 



Po + 



A 



ttA (x — vty + 6^ 
X — vt 



± (A - 1) tan-i 



6 r] 



1 TTlf^t 

2 /i 



vr/iApo 



771 



1 + 



A-1 



1/2 



(73) 



Thus the velocity and width are related to each other. From Eq. (^) we see 
that |f I > Vs and any value of 6po is allowed for A > 1. For A < 1, we must 
have 6po > (1 — A)/7rA; then |f | < Vg- These ranges of velocity agree with the 
exact results known for the particle and hole respectively, as discussed after 
Eq. (^. The identification with particle for A > 1 and hole for A < 1 may 
be justified by considering the sign of the integrated density for the soliton, 

A-1 



dx [ p{x, t) - Po ] 



A 



(74) 



The magnitude of this number is generally not an integer; the physical mean- 
ing of this is not clear to us. It is interesting to note at this point that if 
we perform the scaling (^) which eliminates A, then the redefined soliton 
number is 1 for all A. Hence (|73|) describes a one particle solution in the 
redefined theory. 

The momentum and energy (obtained after subtracting the background 
value) of the soliton ( |75D are given by Eqs. (|TBp and (^). 

A-1 



P 

E 



A 
A-1 



mv 



A ' 2 
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mv 



p). 



(75) 



This dispersion relation does not agree with the exact dispersion relations 

given in (^ and (^. It therefore seems that the interpretation of solitons as 

particles or holes has some difficulties which need to be resolved. 

It is interesting to observe that we can also go in the opposite direction 

and recover the large amplitude waves by superposing a number of single 

soliton solutions in a periodic way |T^ . For this purpose, it is useful to recall 

that 

2^ -y I ^ sinhg 

n£^oo (2vrn + kx)"^ + a^ cosh a — cos {kx) 
It is worth noting that for A > 1, there is no lower bound on po! in 

particular, we can set po = 0. We then get a new solution corresponding to a 

stationary and isolated soliton with no background density. This may also be 

seen for the static case from Eq. (|29|) with no external harmonic potential, 

Ap„ + A^*=0. (77) 

The solution of this is found to be 

/ N A — 1 h ,„„, 

ttA x^ + 0'^ 



with eigenvalue Eq = 0. We may now boost this solution using (pi]) . The 

general solution is therefore 

. , A-1 b 

p{x,t) = 



vrA {x — vty + IP' 



mv 



h 



1 



e{x,t) = -^ X vt , (79) 



2 



where the width and velocity are now independent parameters. The particle 
number, momentum and energy of this soliton are given by 
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p = 


A-1 
A "'^' 


R — 


A — 1 mf ^ 



A 2 - («•" 



Note that if we had used the scaled variables given by Eq. ( P5D , the normal- 
ized dx J p{x) = 1, with the momentum and energy like a classical particle. 
Even though it is an exact solution of the collective field theory, it not a 
genuine solution of the CSM. This is because the collective field theory is 
meaningful only for large A^. 

Finally, we note that we may rewrite Eq. (^) for the scaled density p as 

Formally, this equation has the same form as the steady-state Coulomb gas 
model of Dyson [0. In the diffusion problem, it is known as the Smolu- 
chowski equation with a singular kernel, and describes the Brownian motion 
of a particle immersed in a fluid, with friction-limited velocity. A description 
of this equation is given by Andersen and Oppenheim [^. The analogous 
single-soliton solution ([78| ) of the equation for the diffusion problem was ob- 
tained by Satsuma and Mimura [Q. These authors also found the soliton 
with the hyperbolic kernel, and the periodic solution appropriate for the 
Sutherland Hamiltonian on a circle |0]. 

7. DISCUSSION 

We have seen that collective field theory is a powerful technique from 
which many properties of the CSM can be derived without having to solve 
the A^-particle Schrodinger equation (|^). We can consider other applications 

26 



of collective field theory. For instance, it should be possible to solve for 
the low-energy excitations of the CSM in a slowly time-varying harmonic 
potential. Analogous calculations have been performed for a trapped Bose- 
Einstein condensate in the three-dimensional problem p4| . 

It is evident that there are several issues which are either not clear or 
beyond the reach of collective field theory. We list some of these problems 
below; they are of course related to each other. 

(a) The dispersion relation of small amplitude waves whose wavelengths are 
comparable to the average particle spacing remains unknown. The difficul- 
ties mentioned in Section 4 seem to suggest that the collective field theory 
discussed in this paper is not complete; perhaps one needs higher derivative 
terms in the energy functional (^) to obtain better results at short wave- 
lengths. 

(b) The interpretation of the large amplitude waves and soliton solutions 
is not clear. The ideal way to resolve these difficulties would be to set up 
a precise correspondence between the collective excitations and the known 



solutions of the Schrodinger equation. Refs. [^, |T^ make some suggestions 
in this direction, but a quantitative mapping is still missing. 

(c) It is not clear why we only get exact soliton solutions which correspond to 
particles for A > 1 and holes for A < 1. It would be desirable to complete the 
story by finding, perhaps numerically, solutions corresponding to particles 
for A < 1 and holes for A > 1. 

(d) Finally, it would be very interesting to quantize the collective field theory 
and study it more carefully than we have done in Section 4. This may 
lead to an alternative way of deriving the oscillating terms in the dynamical 
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correlation functions. 
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